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I. INTRODUCTION 

The invariant masses of the stable charmed hadrons 
are currently reported by the Particle Data Group (PDG) 
with a precision of about 0.5-0.6 MeV/c 2 [1]. The best in- 
dividual measurements have a statistical and systematic 
precision of about 0.5McV/c 2 and use data samples of 
a few hundred events. The BABAR data contains large 
samples of many charmed-hadron decays and, due to 
the excellent momentum and vertex resolution in BABAR, 
many of the decay modes can be reconstructed with an 
event- by-event mass uncertainty of a few MeV/c 2 . We 
can therefore significantly improve the precision on the 
charm-hadron mass measurements. 

In this analysis we present a precision measurement 
of the Af mass. The measurement is based on the re- 
construction of the decay modes Af — > AK®K + and 
A+ -> E°K°K+ [2]. Because almost all of the A+ in- 
variant mass in these decays results from the well-known 
rest-mass values of the A+ decay products, the systematic 
uncertainty in the reconstructed mass is significantly re- 
duced compared to the precision obtained in other decay 
modes. Large samples of A+ — > pK~n + and Af — > pK s 
decays are used for cross-checks and studies of systematic 
uncertainties. 



II. ANALYSIS METHOD 

The measurement of invariant mass relies on precise 
and unbiased measurements of particle three-momenta. 
If the chosen decay mode contains any photons, unbiased 
energy and position measurements in the electromagnetic 
calorimeter (EMC) are also necessary. The momentum 
measurement depends on a well-aligned detector, precise 
knowledge of the magnetic field and material distribu- 
tion in the tracking volume, and mass-dependent cor- 
rections for the energy loss of charged particles passing 
through the detector material. All three requirements 
are met to a large degree by the BABAR detector and its 
event reconstruction algorithms. We quantify residual 
systematic effects in studies of various control samples. 
To minimize such contributions to the systematic uncer- 
tainty, we choose to measure the A+ mass by using decay 
modes with a low Q-value, where the Q-value for a decay 
a — > b + c + . . . is defined as 

Q = m(a) — m(b) — m(c) — . . . 

Uncertainties related to track reconstruction, such as 
those involving energy-loss correction or magnetic field 
strength, tend to scale with the Q-value. 

The main signal mode used in this analysis is A+ — > 
AK°K + , which has a Q- value of 177.9 MeV/c 2 calculated 
using the PDG A+ mass. The branching ratio T(A+ — > 
AK°K+)/T(A+ -» P K-tt+) was measured by CLEO to 
be 0.12±0.02±0.02 [3]. BELLE observed that about 25% 
of these decays proceed through E(1690)°K + [4]. Since 



the £'(1690) baryon has a width of several MeV/c 2 and 
its mass is not well-known, it does not help in constrain- 
ing the A+ mass. We reconstruct the A and K® only in 
their charged decay modes, A — ► pn~ and — ► n + n~ , 
which account for about 44% of A+ — ► AK^K + decays. 
To obtain the most precise A+ mass value, the long-lived 
particles are reconstructed with their mass values con- 
strained to their respective PDG values. 

A second low Q-value A+ decay mode, A+ — > 
£°K®K + , has also been studied. With the Q- value of 
100.9 MeV/c 2 it has the potential for an even lower sys- 
tematic uncertainty. The results of this study indicates 
that the branching fraction for this decay mode is signif- 
icantly smaller than for the A+ — > AK° S K + decay mode, 
which makes the statistical precision worse than the total 
precision of the 4+ — > AK®K + mode. The decay also re- 
quires the measurement of a photon from the S° — > A7 
decay and this introduces an additional systematic un- 
certainty. 

For A+ — > AKgK + decays, the event-by-event uncer- 
tainty on the invariant mass is about 2MeV/c 2 ; hence, 
with a few thousand reconstructed signal events, the sta- 
tistical uncertainty is below 50 keV/c 2 and the dominant 
uncertainty is systematic. The major part of the sys- 
tematic uncertainty is estimated directly from the data 
by redoing the track fits with different assumptions on 
the amounts of detector material or the magnetic field 
strength, and measuring how much the invariant mass 
value changes with each assumption. The candidate se- 
lection criteria have been optimized on simulated events 
to minimize the expected systematic uncertainty. 

To check that the procedure for estimating the system- 
atic uncertainty is reasonable, we study large samples of 
A and K® decays. The goal is to ensure that the mea- 
sured mass values are consistent with the PDG values 
within the systematic uncertainty we estimate using the 
same procedure as for the signal mode. Large samples 
of 4+ — > pK~n + and A+ — > pK s decays are also used 
and give invariant mass measurements that are consistent 
with the signal when their larger systematic uncertainties 
are taken into account. 



III. THE BABAR DETECTOR 

The BABAR detector is described in detail elsewhere [5] . 
The momenta of charged particles are measured with a 
combination of a five-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) in a 1.5-T solenoidal 
magnetic field. The momentum resolution is measured to 
be a(p T )/p T = 0.0013(p T / GeV/c)®0.0045. A detector of 
internally reflected Cherenkov radiation (DIRC) is used 
in charged particle identification. Kaons and protons are 
identified with likelihood ratios calculated from dE/dx 
measurements in the SVT and DCH, and from the ob- 
served pattern of Cherenkov light in the DIRC. A finely 
segmented CsI(Tl) electromagnetic calorimeter (EMC) is 
used to detect and measure photons and neutral hadrons, 



7 



and to identify electrons. The instrumented flux return 
(1FR) contains resistive plate chambers for muon and 
long-lived ncutral-hadron identification. For event sim- 
ulation we use the Monte Carlo generator EVTGEN [6] 
with a full detector simulation that uses GEANT4 [7] . 

The most critical component of this analysis is the 
quality of the track reconstruction. In order to max- 
imize tracking efficiency, the track-finding algorithm is 
based on tracks found by the trigger system and by stand- 
alone track reconstruction in the SVT and in the DCH. 
Once a track has been found, the track parameters are 
determined using a Kalman filter algorithm [8], which 
makes optimal use of the hit information and corrects 
for energy loss and multiple scattering in the material 
traversed and for inhomogencitics in the magnetic field. 
The material-traversal corrections change the track mo- 
mentum according to the expected average energy loss 
and increase the covariance for track parameters to ac- 
count for both multiple scattering and the variance in 
the energy loss. The energy loss depends on the particle 
velocity, therefore each track fit is performed separately 
for five particle hypotheses: electron, muon, pion, kaon, 
and proton. A simplified model of the BABAR detector 
material distribution is used in the Kalman filter algo- 
rithm in order to maintain reasonable execution time. 
The main layers of material traversed by a particle origi- 
nating from the interaction point are the beam pipe at a 
radius of 2.5 cm, consisting of about 1.4-mm of beryllium 
and 1.5-mm of cooling water; five layers of 300-/iin-thick 
silicon detectors at radii from 3.3 cm to 15 cm; a 2-mm- 
thick carbon-fiber tube at a radius of 22 cm supporting 
the SVT and beam-line magnets; and the inner wall of 
the DCH at a radius of 24 cm, which is a 1-mm thick 
beryllium tube. Part of the support structure for the 
silicon detectors is modeled by increasing the thickness 
of each layer in the silicon detector by 60 /im, while the 
effect of the DCH gas is modeled as a series of discrete 
material contributions. Detailed knowledge of the mag- 
netic field is also essential to the track reconstruction. 
This is discussed in Section VII B below. 



V. EVENT SELECTION 

A. A+ -» AK^K+ Selection 

The A+ — > AK®K + signal is reconstructed using only 
the charged two-body decay modes of the A and K° s 
hadrons. We form A candidates from two tracks, one 
of which must be identified as a proton, and, after fitting 
to a common vertex, we require the combined invariant 
mass to be between 1106 and 1125 McV/c 2 . For surviv- 
ing candidates, the two tracks are fit to a common vertex 
with the invariant mass constrained to the PDG A mass. 
The probability of this mass-constrained vertex fit is re- 
quired to be above 10 -3 . Similarly, a candidate is 
formed from two tracks, neither of which belongs to the 
A candidate, with a combined invariant mass between 460 
and 530MeV/c 2 . For surviving K® candidates, the two 
daughter tracks are fit to a common vertex with the mass 
constrained to the PDG K® mass and the fit probability 
is required to be higher than 10~ 3 . The A and K® can- 
didates are then combined with a fifth track, identified 
as a charged kaon, in a fit to a common vertex to form 
a At candidate. The At candidate must have an invari- 
ant mass between 2250MeV/c 2 and 2330MeV/c 2 . The 
probability of the vertex fit should be greater than 10~ 3 . 
To suppress combinatorial background, the signed decay 
length of a K® candidate is required to be larger than 
three times its estimated uncertainty. The signed decay 
length is defined as the distance between the At and K® 
candidate along the K® momentum in the CM frame. To 
further suppress background, which results mainly from 
B decays, the momentum (p*) of the A£ candidate in 
the e + e~ center-of-mass frame (CM) is required to be at 
least 2 GeV/ c. This requirement also helps to reduce sys- 
tematic uncertainties that affect mainly low-momentum 
tracks. The selection efficiency, not including branching 
fractions, is about 15% for A+ — > AK° S K + decays with 
CM momentum larger than 2 GeV/ c. The background 
is suppressed sufficiently to not be an issue for the At 
mass measurement. 



IV. DATA SAMPLE 

The data sample used for the At mass measurement 
comprises an integrated luminosity of 232fb _1 collected 
from e + e~ collisions at or 40 MeV below the T(4S) res- 
onance. For the studies of A and decays only a small 
subsample of the data is used due to the high produc- 
tion rates of these hadrons. Studies of simulated events 
are performed using Monte Carlo samples of generic 
e + e~ — > BB and e + e~ — > qq(q = u,d,s,c) continuum 
events with an integrated luminosity equivalent to 240 to 
275 fb" 1 . More than 230000 simulated /1+ -> AK° S K+ 
decays and 60000 simulated At — > S°K^K + decays are 
used for studies of systematic uncertainty. 



B. At -> S°K S K+ Selection 

The At — > £°K®K + mode is reconstructed from 
U° — > Aj, and A and K® hadrons decaying into two 
charged particles. The A and candidates are formed 
in the same way as in the At -> AK°K + selection. A 
S° candidate is formed by combining a A candidate with 
a photon and requiring the combined invariant mass to 
be between 1184 and 1196MeV/c 2 . A photon candidate 
is defined as an energy deposit in the EMC of at least 
30 MeV that is not associated with any track and has 
a lateral moment [9] of its shower energy deposition of 
less than 0.8. If several photons can be combined with 
a A to form U° candidates, only the candidate with the 
most energetic photon is retained in order to avoid double 
counting. The S° candidates are fit with their mass con- 
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strained to the PDG mass and are combined with K ° and 
K + candidates to form yl+ candidates that must satisfy 
the same invariant mass, vertex probability, and p* re- 
quirements as j4+ — > AK° S K + candidates. The selection 
efficiency, not including branching fractions, is estimated 
to be about 8% for A+ -> S°K°K+ decays with A+ CM 
momentum larger than 2 GeV/c. 

From simulation we expect that 18% of the — > 
S°K^K + decays are reconstructed with the wrong pho- 
ton in the S° candidate. Due to the low energy of the 
photon and the mass constraint on the S° candidate, 
these A+ candidates still have the correct mass on aver- 
age, but the mass resolution of these candidates is signif- 
icantly degraded. 
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FIG. 1: Invariant mass distribution for ylj — » AK%K + can- 
didates. The lower part of the figure shows the normalized fit 
residuals. The dashed line indicates the present PDG value. 



j4+ — » pA"~7r + candidates are formed from three tracks 
identified as a proton, a kaon, and a pion with a com- 
bined invariant mass between 2240 and 2330 MeV/c 2 . 
The tracks are fit to a common vertex and the probabil- 
ity of the vertex fit is required to be greater than 10~ 3 . 
The signal-selection efficiency is about 42% and depends 
on the At momentum. 



D. A+ -* pK% Selection 

For the yl+ — » pK° s decay mode, the K® candidates 
are required to satisfy the same criteria as in the A£ — » 
AK®K + mode, but with the further constraint that the 
decay angle 9 of the ir + in the K® rest frame with respect 
to the K° s line-of-flight must satisfy | cos 6»| < 0.97. This 
removes contamination from 7 conversions. The K® can- 
didates are combined with tracks identified as protons in 
a fit to a common vertex and the resulting Aj candidates 
are required to have an invariant mass between 2240 and 
2330 MeV/ c 2 . The probability of the vertex fit is required 
to be above 10 -3 . The signal selection efficiency is about 
41% and depends on the A+ momentum. 



E. A — > pn~ and K° 



Selection 



For the cross check studies in Section VII, large sam- 
ples of A and K® decays are reconstructed using similar 
criteria as those used to select A and K® candidates for 
A£ — ► AKgK + decays. The mass-constrained vertex fits 
are replaced with geometric vertex fits with the require- 
ment of a fit probability greater than 10 -2 . For both 
A and candidates the signed decay length calculated 
with respect to the e + e~ interaction point is required to 
be larger than three times its uncertainty. 



TABLE I: Fit parameter values for the Aj -» AK° S K + and 
At — ► S°K°K + signals together with the half- width at half- 
maximum (HWHM) calculated from these values. A correc- 
tion to the mass values for underestimated energy loss has not 
been applied. 



Parameter 

Fitted mass ( MeV/c 2 ) 
Signal yield (events) 
Narrow width (MeV/c 2 ) 
Broad width (MeV/c 2 ) 
Narrow fraction 
HWHM (MeV/c 2 ) 



At 



AK° q K+ At -> S°K° q K+ 



2286.44±0.04 
4627±84 
2.08±0.07 
6.39±1.22 
(83±3)% 
2.55±0.06 



2286.29±0.18 
264±20 
2.04±0.18 



2.41±0.22 



VI. SIGNAL FITS 

The invariant mass distribution for the A+ — > AK®K + 
candidates is shown in Fig. 1. A clear A* signal peak is 
observed. A binned maximum likelihood fit to the mass 
distribution is performed using a sum of two Gaussians 
with a common mean for the A+ signal function. The 
background is described by a linear function as suggested 
by simulation studies. The fit parameter values are given 
in Table I. Note that the uncertainty on the mean mass 
is statistical only and a correction for underestimated 
energy loss described in the next section has not been 
applied to the fitted mass. 

The invariant mass distribution for the A~£ — > 
£°K®K + candidates is shown in Fig. 2. A small but 
significant Aj signal peak is observed. The figure also 
shows the expected background under the Aj peak from 
A+ — » S°KgK + decays with a correct A but a wrong 
photon used in candidate selection. A binned maximum 
likelihood fit of the mass distribution is performed using a 
single Gaussian for the A+ signal. The background is de- 
scribed by a linear function. The wrongly-reconstructed 
A+ — > E°KgK + candidates are absorbed into the signal 
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FIG. 2: Invariant mass of reconstructed A+ -> E°K° S K + 
candidates. The hashed histogram is the expected contribu- 
tion from signal events with a wrongly assigned photon. The 
lower part of the figure shows the normalized fit residuals. 
The dashed line indicates the present PDG value. 
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FIG. 3: Invariant mass distribution of candidates used in 
cross-checks of systematic uncertainty, a) A — > pn~ de- 
cays, b) Kg — > tt + tv~ decays, c) At — > pK~ir + decays, 
d) At — > pKs decays. The dashed lines show the PDG mass 
values for these hadrons. A correction to the mass distribu- 
tions for underestimated energy loss has not been applied. 



and background because they peak at the A+ mass. The 
fit parameter values are given in Table I. 

The invariant-mass distributions for the four control 
modes are shown in Fig. 3. The A and K® signals are fit 
to a sum of three Gaussians with common mean, while 
the two A+ signals are fit to a sum of two Gaussians 
with common mean. The background in all four cases is 
modeled with a second-degree polynomial. The fit yields, 
mass, and signal RMS values are listed in Table II. The 
fitted mass values for the A and Kg are significantly be- 
low the PDG values. This is mainly due to an underes- 
timation of the energy loss in the detector material and 
is described in more detail in the next section. 



VII. CROSS-CHECKS 

The four control mode samples discussed above are 
used to understand the accuracy to which particle masses 
can be measured in BABAR. 



A. Material Dependence 

The main systematic uncertainty on the A+ mass 
comes from uncertainties in the energy-loss correction in 
charged particle tracking. The low A and K® fitted mass 
values given in Table II indicate that the energy-loss cor- 
rection may be underestimated. The long lifetimes of A 
and K° s hadrons enable us to study this in more detail. 
Figure 4 shows the A and fitted mass values as a func- 
tion of the radial distance from the interaction point to 
their decay point. The further the decay point is from the 
interaction point, the less material the charged daughter 
particles traverse, and so energy-loss corrections become 
less significant. The deviation from the PDG value is seen 
to be largest for decays closest to the interaction point, 
thereby strongly indicating that the underestimation of 
the mass values is due to insufficient material corrections. 
The lower fitted mass values at radii of 12 cm, 21 cm, and 
23 cm coincide with vertices reconstructed inside or very 
near material sites. The same effect is observed in the 
reconstruction of simulated data. This is related to de- 
tails of the incorporation of energy-loss corrections into 
the track-fit procedure in such circumstances. 

The effect of increasing the amount of material as- 
sumed during the track reconstruction has been studied 
using several different scenarios. Figure 4 shows what 
happens if the material density in all parts of the SVT is 
increased uniformly by 20%. This is a gross simplifica- 
tion; however it removes most of the dependence on decay 
radius and moves the fitted hadron masses closer to their 
PDG values. The K® mass is, however, still consistently 
low compared to the PDG value by about O.I5MeV/c 2 . 

Another way to investigate the energy-loss correction 
is to study the fitted mass as a function of particle mo- 
mentum. The lower the momentum of a charged particle, 
the more significant the energy-loss corrections become. 
Candidates are reconstructed from multiple tracks, but 
to simplify this investigation, the mass value is studied 
as a function of net candidate momentum calculated in 
the laboratory frame. Figure 5 shows the fitted mass 
value for A~£ — ► pK~ir + and A+ — ► pK® candidates as 
a function of the laboratory momentum. Above 3 GeV/c 
the fitted mass value reaches a limit that is significantly 
above the PDG, while below 3 GeV/c it falls by more than 
800 keV/c 2 as the momentum decreases. Similar behavior 
is observed for A and K® hadrons, but the limiting value 
is reached at about 2 GeV/ c. The decrease in mass is only 
about 60keV/c 2 for A baryons and about 500keV/c 2 for 
Kg mesons. 

Increasing the assumed material density in the SVT 
by 20% is seen to greatly reduce the radial mass depen- 
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TABLE II: Fit parameter values and calculated HWHM for the four decay modes used for cross-checks. A correction to each 
mass value for underestimated energy loss has not been applied. 



Parameter 


A — > pir 


K° s - n+n- 


A+ -^pK-n+ 


A+ -> P K° S 


PDG mass (MeV/c 2 ) 


1115.683±0.006 


497.648±0.022 


2284.900±0.600 


2284.900±0.600 


Fitted mass (MeV/c 2 ) 


1115.660±0.001 


497.305±0.002 


2285. 845±0. 013 


2285.876±0.023 


Signal yield (events) 
HWHM (MeV/c 2 ) 


3192700±5800 


2463900±4900 


1449300±5300 


243700±1900 


0.853±0.002 


2.715±0.005 


5.147±0.014 


5.613±0.046 




FIG. 4: Fitted mass values for A — ► pir~ (left) and K s — > -k + tt~ (right) candidates minus the PDG value as a function of decay 
vertex radius. The solid circles are for track reconstruction with the normal amount of detector material. The open squares 
are for track reconstruction assuming 20% more material in the SVT. The dashed lines show the ±lcr uncertainty on the PDG 
value for the masses. 



dence in the A and K® control samples and the momen- 
tum dependence in all control samples. We therefore 
apply this change to the reconstruction of the two low 
Q-value modes in order to obtain a more accurate mass 
measurement. However, since the K® mass still is not in 
agreement with the PDG value, we will use the largest 
observed variation in the fitted mass when we vary the 
material model as an estimate of the systematic uncer- 
tainty. 



B. Magnetic Field Dependence 

The momentum measurements, and thus the mass 
measurements, depend critically on the magnetic held. 
The main component of the magnetic field is the solenoid 
field, which has an average value of 1.5 T parallel to the 
beam axis. The field was mapped very precisely with 
movable Hall probes before the detector was installed; an 
NMR probe measured the absolute field strength. The 
field strength in the tracking volume is estimated to be 
known to an accuracy of about 0.2 mT. The second most 
significant field component comes from the permanent 
magnets used for the final focusing and bending of the 
beams. The magnet closest to the interaction point is 
about 20 cm away. The fringe fields in the tracking vol- 
ume from the magnets are weak and have been measured. 



More uncertain is the contribution to the magnetization 
of the permanent magnet material due to the solenoid 
field. This effect is measured only at a few specific points 
with Hall and NMR probes and we therefore rely on a 
finite element calculation to estimate the magnetization 
effect elsewhere. This model depends on the permeabil- 
ity of the magnets. These magnets are made of a SmCo 
alloy, which has a measured permeability /i = 1.07 in the 
direction transverse to the solcnoidal field. However this 
is an average over many samples, which range from 1.04 
to 1.10. Furthermore in the direction of the solenoidal 
field the permeability of the SmCo elements is measured 
to be about 5% larger. Therefore there is a significant 
uncertainty on this component of the field. The longitu- 
dinal field component at the interaction point from these 
magnets is about +9 mT. This increases the curvature of 
charged particle trajectories. 

We vary the assumed solenoid field strength by 0.02% 
for the systematic uncertainty studies. The magnetiza- 
tion field is varied by 20% in order to account for differ- 
ences between the direct field measurements and the per- 
meability measurements. Fig. 6 shows the effect on the 
fitted At — » pK~ir + mass value as a function of momen- 
tum. The fitted mass value is seen to be shifted by the 
same amount independent of the momentum. The same 
is true for the other control samples, in particular for the 
A and K® mass dependence on decay vertex radius. The 
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FIG. 5: Fitted mass value for At — > pK~iv + (left) and At — ► pK° (right) candidates minus the PDG value as a function of At 
laboratory momentum. The solid circles are for track reconstruction with the normal amount of detector material. The open 
squares are for track reconstruction assuming 20% more material in the SVT. The dashed line shows the +la uncertainty on 
the PDG value of the At mass. 



,3-2000 




5 6 
P, ab (GeV/c) 



FIG. 6: Fitted mass of At — * pK~ir + candidates minus the 
PDG value as a function of A* laboratory momentum. The 
additional energy-loss correction has not been applied. The 
solid circles are for track reconstruction with the normal mag- 
netic field. The open squares are for track reconstruction as- 
suming a 0.02% higher solenoid field. The open triangles are 
for track reconstruction assuming 20% higher magnetization. 
The dashed line shows the +la uncertainty on the PDG value 
of the At mass. 



control samples therefore provide little guidance on the 
magnetic field uncertainty. 



C. <(> Dependence 

Studies of the control samples reveal a significant de- 
pendence of the fitted mass value on the azimuthal angle 
<j> of the hadron candidate momentum at the origin. This 
effect is not seen in simulated events. The <f> dependence 
is shown in Figure 7 for the K® and Af — > pK~n + sam- 



ples. The effect is seen to be roughly anti-symmetric in 
</> with the lowest fitted mass at cfi = §, corresponding 
to upward going hadrons, and the largest fitted mass at 
<f>=*Z-, corresponding to downward going hadrons. 

We have not been able to identify the source of the <j> 
dependence, but we can estimate the potential impact on 
the At mass measurement. The magnitude of the <j> vari- 
ation increases with the momentum of the reconstructed 
hadron. This shows that the variation is not due to some 
asymmetry in the material distribution that was unac- 
counted for in reconstruction, which would influence low 
momentum particles the most. The <f> dependence is also 
observed when only the DCH is used for reconstructing 
the tracks. This indicates that the source may be related 
to the internal alignment of the DCH, which is based on 
end-plate surveys done during construction of the cham- 
ber. 

The <f) dependence in the fitted mass can be repro- 
duced qualitatively in simulated events by introducing an 
explicit 4> dependence in the measured track momenta. 
To get the largest <j> dependence at high momentum, 
the change has to be introduced as a bias in the track 
sagitta measurements. To a good approximation this cor- 
responds to changing the transverse track momentum pr 
according to 



1 1 r • i 

= h o sin <p. 

Pt' Pt 



(1) 



The constant S is chosen to reproduce the observed mag- 
nitude of the <j> dependence. The chosen value corre- 
sponds to a change in track sagitta of 11 /um, where the 
full sagitta for a track with 1 GeV/ c transverse momen- 
tum is about 3 cm. Applying this modification to the 
A+ -» AK° S K+ and A+ -» S°K°K+ Monte Carlo sam- 
ples introduces a <j> dependence in the fitted mass values 
with an amplitude of 55keV/c 2 and 30keV/c 2 , respec- 
tively. When averaging over all 4>, the fitted mass value 
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FIG. 7: Fitted mass value of K° — > -k + tt~ (left) and A~£ — > pK~ir + (right) candidates minus the PDG value as a function of the 
azimuthal direction. The additional energy- loss correction has not been applied. The dashed lines show the ±l<r uncertainty 
on the PDG value of the masses. 



does not change, but given that the source of the de- 
pendence is not understood, we use these amplitudes as 
estimates of systematic uncertainty. 

VIII. SYSTEMATIC UNCERTAINTIES 

The major sources of systematic uncertainty (energy- 
loss correction, magnetic field, and cf> dependence) have 
been described in detail in the previous section. All of 
the known systematic uncertainties are listed in Table III. 
For the systematic uncertainty on the energy-loss correc- 
tion we use the observed change in mass when the mate- 
rial density is increased uniformly by 10% in the track- 
ing volume, which is a slightly larger change than the 
one from increasing the density in the SVT by 20%. The 
fit-procedure uncertainty summarizes the variation in the 
fitted mass value when the shapes of the signal, the back- 
ground, and the choice of binning are varied. Possible bi- 
ases from the internal alignment of the SVT are studied 
by applying small distortions to the SVT alignment in 
simulated events. The magnitude of the distortions cor- 
respond to the changes observed in the internal alignment 
between different run periods. For the — > E K®K + 
mode, we vary the EMC energy scale by 5%, but due to 
the low energy of the photon and the mass constraint on 
the 17°, this has little effect on the fitted mass. Finally 
the A, K + and K® masses have uncertainties of 6, 16 and 
22 kcV/ c 2 and the effect of these uncertainties on the A+ 
mass has been estimated. 



A. Energy-Loss Correction 

The fitted A~£ mass values from Section VI need to be 
corrected for the underestimated energy loss. The cor- 
rection is calculated by increasing the material density 
of the SVT by 20%, which is seen in the control samples 



TABLE III: Systematic uncertainty contributions to the A+ 
mass measurements from the A£ — > AK° S K + and A£ — ► 
S°K° S K + samples (in keV/c 2 ). 



At- 


-> AK° S K+ 


A+ -> E°K° S K + 


Solenoid field 


±60 


±30 


Magnetization 


±68 


±29 


Energy-loss correction 


±83 


±50 


4> dependence 


±55 


±30 


Fit procedure 


±38 


±71 


SVT alignment 


±23 


±13 


EMC energy scale 




±46 


Particle masses 


±27 


±58 


Total systematic 


±144 


±126 



to remove most of the momentum and decay-radius de- 
pendence. The corrections are 61 keV/c 2 and 18 keV/c 2 
for the A+ AK° S K+ and A+ -> Z°K°K+ samples, 
respectively. This gives the following results for the A+ 
mass: 

m(A+) AK o K + = 

2286.501 ± 0M2(stat.) ± QAAA{syst.) MeV/c 2 , 
m(A+) E o K o K+ = 

2286.303 ± 0A8l(stat.) ± 0A26(syst.) MeV/c 2 . 

B. Combined Result 

The systematic uncertainties on the two measurements 
are highly, but not fully, correlated. We combine the two 
mass measurements using the BLUE (Best Linear Unbi- 
ased Estimate) technique [10]. Besides the statistical un- 
certainty, we consider the fit procedure uncertainty and 
uncertainties related only to the A+ — ► £°K®K + mode 
to be uncorrelated, while the remaining systematic uncer- 
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TABLE IV: Fitted At mass values and major contributions to TABLE V: Fitted A and K$ masses and the major contribu- 
the systematic uncertainty for the At — > pK~TT + and At — > tions to the systematic uncertainty for candidates from the 
pK® samples used for cross-checks (in MeV/c 2 ). At — * AKgK + sample (in MeV/c 2 ). 





A+- 


-> pK 7T + 


A+ - pK° s 




A 


K° 


Fitted Mass 


2286. 182±0. 018 


2286.216±0.034 


Fitted Mass 


1115.657±0.014 


497.359±0.040 


Corrected Mass 


2286.393±0.018 


2286.361±0.034 


Corrected Mass 


1115.679 ±0.014 


497.560±0.040 


Solenoid field 




±0.181 


±0.196 


Solenoid field 


±0.009 


±0.068 


Magnetization 




±0.207 


±0.222 


Magnetization 


±0.005 


±0.054 


Energy-loss correction 




±0.278 


±0.199 


Energy-loss correction 


±0.040 


±0.242 


cj> dependence 




±0.217 


±0.236 


4> dependence 


±0.006 


±0.059 


Total systematics 




±0.447 


±0.428 


Total systematic 


±0.041 


±0.264 



taintics arc 100% correlated. The correlation coefficient 
for the two measurements is estimated to be 0.355. The 
combined mass result is 

m(A+) = 2286.46 ± 0.14 MeV/c 2 . 



C. Mass Cross-checks 

From the two large-Q-value Af data samples, we ob- 
tain measurements of the A£ mass that can be com- 
pared to our more precise measurements from the At — > 
AK°K+ and A+ S°K° S K+ samples. To keep the 
systematic uncertainty low, we use only Af — > pK~n + 
and At — ► pK® candidates with laboratory momentum 
above 3 GeV/c, as Figure 5 shows that those candidates 
have less dependence on the assumed amount of detector 
material. The resulting mass value for each decay mode 
from fitting the invariant-mass spectra with the sum of 
two Gaussian distributions with common mean is given 
in Table IV. The central values are corrected for the shift 
in mass observed when the material density of the SVT is 
increased by 20%. These corrections are 211kcV/c 2 and 
145keV/c 2 for the A+ — > pK^n + and A+ — > pK° modes, 
respectively. 

Table IV also lists the major systematic uncertainties. 
The solenoid and magnetization fields are varied as for 
the low-Q-value modes. For the energy-loss correction, 
we compare the effect of increasing the density of the 
SVT by 20% to the effect of increasing the density of ma- 
terial in the tracking volume by 10%, taking the larger 
change from the standard reconstruction as an estimate 
of the uncertainty. For both decays the larger effect is the 
10% material change in the full tracking volume. The <fi 
dependence is estimated by introducing a 4> dependence 
in the simulation as described by Eq. (1). The uncer- 
tainty is the maximum change in mass introduced. The 
results 

m( A t)pK-ir+ = 

2286.39 ± 0.02(stat.) ± 0A5(syst.) MeV/c 2 , 
™{A+) pK o = 

2286.36 ± 0.03(stat.) ± 0A3(syst.) MeV/c 2 . 



arc in good agreement with our main result, but have 
larger systematic uncertainties. 

The At — > AKgK + sample is sufficiently large that 
we can fit the A and K® mass distributions for candi- 
dates that are combined to form the Af — ► AK®K + . In 
order to fit the A and mass distributions, the mass 
constraint is removed from the candidate reconstruction 
and the resulting invariant mass spectra are fit with a 
double- Gaussian signal shape and a linear background. 
The resulting mass values and systematic uncertainties 
are listed in Table V. The mass correction for the under- 
estimated energy loss and the systematic uncertainties 
are estimated as for the A+ — > pK~n + and A+ — ► pK° 
samples. The final results after the energy-loss correc- 
tion, 

m(A) = 1115.68 ± 0.01(stat.) ± 0.04(sysi.) MeV/c 2 , 
m{K° s ) = 497.56 ±0.0A(stat.) ±0.26{syst.)MeV/c 2 , 

are in agreement with the PDG 

m PDG {A) = 1115.683 ± 0.006 MeV/c 2 , 
m PDG {K° s ) = 497.648 ± 0.022 MeV/c 2 . 

Since the A and K® candidates are the same candidates 
used in the final At sample, the agreement with the PDG 
mass values gives further confidence in the Af mass re- 
sult. 



IX. SUMMARY 

We have presented a precision measurement of the At 
mass using the low-Q-value decay modes Af — > AK®K + 
and At — > £°K®K + in order to minimize systematic 
uncertainty. The measured mass in the two modes is 

m(A+) AK o K+ = 

2286.501 ± 0M2(stat.) ± 0AU(syst.) MeV/c 2 , 

2286.303 ± 0A8l(stat.) ± 0A26{syst.) MeV/c 2 . 

Combining these measurements, taking the correlated 
systematics into account, the final result for the At mass 
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is 

m(A+) = 2286.46 ± 0.14MeV/c 2 . 

This result is in agreement with the mass values mea- 
sured in other A£ decay modes, although those are sub- 
ject to large systematic uncertainty. The systematic un- 
certainty has been cross-checked using large data samples 
of A, K® and A£ decays. The studies have shown that 
there is an underestimation of the energy-loss correction, 
and a dependence on azimuthal angle in the standard 
BABAR track reconstruction. The impact on the mass 
measurement has been taken into account in the corre- 
sponding estimates of systematic uncertainty. 

This A£ mass measurement is the most precise mea- 
surement of an open charm hadron mass to date and is an 
improvement in precision by more than a factor of four 
over the current PDG value of 2284.9 ± 0.6 MeV/c 2 . Our 
result is about 2.5a higher than the PDG value, which 
is based on several high Q-value decay modes, mainly 
A+ — ► pK~ir + decays. Theoretical calculations of the A+ 
mass, such as those based upon lattice QCD [11] or ad- 
vanced potential models [12], currently have significantly 
larger uncertainty than the experimental result presented 
here. 
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